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FINAL  REPORT 


Introduction 

The  principal  objectives  of  this  research  program  are  to  investigate  the 
spatial  and  temporal  variability  in  infrared  (IR)  radiance  produced  by  atmospheric 
gravity  waves.  The  investigations  will  be  divided  into  the  following  four  parts: 

(1)  The  spatial  and  temporal  variations  in  concentrations  of  [01,  [NO] , 

[ CO^ 1  and  [I^O]  produced  by  gravity  waves  for  different  exospheric  temperatures 
corresponding  to  cases  A,  C,  F,  G  and  H  in  the  AFGL  report  (AFGL-TR-81- 
0207) . 

(2)  The  spatial  and  temporal  variations  in  the  infrared  volume  emission 
rates  of  some  of  the  above  mentioned  species. 

(3)  The  development  of  a  dispersion  relation  for  analyzing  short  period 
gravity-wave  induced  infrared  variance. 

(4)  The  possible  production  of  short-period  (of  the  order  of  minutes)  gravity- 
wave  oscillations  which  can  be  either  a  primary  wave  or  a  secondary  gravity 
wave  produced  by  primary  gravity  waves  with  large  amplitudes  (20%  or  greater) . 
Such  basically  horizontal  "surface-like"  waves  occur  often  in  the  mesosphere 
with  horizontal  wave  length  in  the  tens  of  kilometers  and  horizontal  phase 
velocities  primarily  in  the  few  tens  of  meters  per  sec.  These  waves  tend 

to  produce  very  significant  spatial  and  temporal  variations  in  infrared 
radiance. 

The  possible  connection  between  oscillations  in  atmospheric  emissions  and 
atmospheric  waves  was  noted  as  far  back  as  1957  by  Krassovsky  (1957).  Later, 
these  effects  were  observed  and  analyzed  by  a  number  of  investigators  [e.g.  Okuda, 
(1962);  Silverman,  (1962);  Barbier,  (1964,  1965);  Weill  and  Christophe-Glaume, 
(1967);  Dachs,  (1968);  Andrews,  (1976);  Dyson  and  Hopgood,  (1978)  for  both  the 
6300A  01  and  the  5200  A  NI.  The  presence  of  gravity  waves  in  OH  infrared 


emissions  has  been  examined  by  Krassovsky  (1972),  Krassovsky  and  Shagaev  (1974,1977), 
Krassovsky  et  al  (1977),  Armstrong  (1975)  Moreel  and  Herse  (1977)  and  Peterson 
(1979).  Theoretical  treatments  for  the  effects  of  gravity  waves  on  airglow  have 
been  made  by  Porter  et  al  (1974)  for  6300A  01  and  5200  A  NI  and  Hatfield  et  al 
(1981)  for  OH  infrared  emissions. 

In  general,  the  treatment  of  gravity-wave  effects  on  airglow  requires  an 

analysis  of  the  effects  of  gravity  waves  on  the  individual  atmospheric  constituents 

which  produce  the  airglow.  Thus,  in  part  (1)  we  analyze  the  effects  of  gravity 

waves  on  the  concentration  of  COj*  NO,  0  and  H^O  using  the  same  basic  technique 

developed  by  Hatfield  et  al  (1981)  for  0^.  For  gravity-wave  models  we  use  the 

models  developed  by  Yu  et  al  (1980).  They  were  developed  for  the  exospheric 

temperatures  T^  equal  to  600°k,  1000°k,  1500°k,  1800°k  and  2200°k,  based  on 

the  1972  COSPAR  atmosphere.  These  gravity-wave  models  are  particularly  suitable 

for  our  purposes,  since  the  exospheric  temperatures  correspond  very  closely  to 

cases  A  (T  =  626. 5°k),  C(T  *  1081. 2°k),  F(T  *  1555. 7°k),  G(T  -  1829. 4°k) 

6x  6x  ex  ex 

and  H(Tgx  »  2176. 8°k)  that  are  being  used  in  the  AFGL  report  (AFGL-TR-81-0207) . 

For  part  (2)  once  the  variations  in  the  concentration  profiles  (produced 
by  gravity  waves)  are  obtained,  the  production  of  infrared  emission  from  these 
species  depend  entirely  on  the  mechanism  which  put  the  species  in  a  particular 
excited  state.  For  OH  infrared  emission,  the  primary  mechanism  is  through  chemical 
reactions,  in  particular  the  reaction  between  H  and  0^  which  produce  the  vibrationall; 
excited  OH  molecules.  For  CO2,  NO,  0  and  HjO,  the  primary  excitation  mechanism 
is  through  collision  (without  reaction)  and,  above  all,  through  absorption  of 
radiation  from  different  sources.  In  all  optical  emission  calculations,  the 
quantity  we  need  to  know  is  the  concentration  of  excited  states.  Once  this  is 
found,  the  volume  emission  rate  can  be  obtained,  usually  by  asstiming  spontaneous 
emission  and  neglecting  stimulated  emissions.  For  low  altitudes,  the  concentration 
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of  excited  states  can  be  found  through  assuming  chemical  equilibrium  (i.e. 
production  rate  of  the  excited  states  is  set  equal  to  the  loss  rate  and  both 
of  which  depend  only  on  reaction  rates)  .  With  this  assumption,  the  concentration 
of  excited  states  can  be  found  algebraically.  This  method  has  been  used  in  the 
AFGL  Report  (AFGL-TR-81-0207) .  At  higher  altitudes,  however,  we  can  no  longer 
assume  chemical  equilibrium,  since  vertical  diffusion  may  become  important.  With 
the  addition  of  this  effect  to  the  loss  rate,  we  end  up  with  a  second  order  self- 
adjoint  differential  equation  for  the  concentration  of  excited  states.  The  volume 
emission  rate  can  be  computed  from  this  concentration  by  using  the  usual  Einstein 
coefficients. 

Of  special  importance  to  infrared  radiance  variance  are  the  short  period 
(of  the  order  of  minutes)  waves  or  "ripples"  which  have  pronounced  spatial  temporal 
variance.  Such  waves  often  (but  not  always)  tend  to  behave  like  surface  waves 
propagating  horizontally  with  relatively  low  phase  velocities  (mostly  less  than 
50  m  sec”*).  These  waves  have  been  detected  at  90  km  by  Peterson  (1979)  through 
their  effects  on  OH  emissions.  To  analyze  these  waves  we  need  a  dispersion  relation 
which  provides  a  formula  for  the  vertical  wave  length,  horizontal  wave  length 
and  the  frequency  of  gravity  waves.  Although  there  are  numerical  dispersion 
relationships  available  (for  reference,  see  Hines  (1974)  and  Francis  (1975)) 
they  are  listed  for  a  few  definite  gravity  wave  periods  and  horizontal  phase 
velocities.  Thus  for  general  use  there  is  only  the  analytic-formula  of  the  Hines 
dispersion  relation  (1960)  which  is  only  valid  for  uniform  isothermal  atmosphere 
with  no  wave  reflection.  We  are  interested,  on  the  other  hand,  in  the  inhomogenieties 
in  spatial  and  temporal  infrared  radiance  produced  by  gravity  waves  in  the 
mesospheric  region  where  the  gravity-wave  amplitudes  grow  to  its  maximum,  while 
the  sharp  temperature  rise  at  the  base  of  the  thermosphere  produces  heavy  wave 
reflections.  Furthermore,  we  are  also  interested  in  short-period  gravity  waves 


with  periods  close  to  the  Brunt  period.  It  is  therefore  necessary  to  generalize 
the  dispersion  relation  to  an  inhomogeneous  atmosphere  which  takes  reflection 
fully  into  account  (see  Papadopoulos  et  al  (1982))  and  remains  valid  as  we  approach 
the  Brunt  frequency.  Such  a  dispersion  relation  will  be  derived  in  this  report 
as  part  of  part  (3). 

For  part  (A)  we  shall  be  primarily  concerned  with  how  a  secondary  short- 


period  oscillation  or  wave  may  be  produced  from  a  primary  long  -period  wave. 

Based  on  the  fact  that  the  velocity  amplitude  of  a  primary  wave  produced  in  the 
lower  atmosphere  must  grow  with  height  until  saturation  of  wave  amplitudes  occur 
in  the  upper  mesospheric  region  where,  at  least  for  the  larger  amplitude  waves, 
the  Richardson's  number  turns  negative  and  instability  occurs. 

Whilst  a  large  variety  of  instability  may  occur  in  the  atmosphere,  it  has 
long  been  assumed  that  an  important  source  of  instability  and  atmospheric  turbulence 
in  the  mesospheric  region  may  be  due  to  the  growth  of  gravity-wave  velocity  amplitudes 
with  height  until  Richardson's  number  falls  below  some  critical  value  (Hodges 
(1967);  Hines  (1974);  Roper  (1979);  Zimmerman  and  Murphy,  (1977);  Manson  and 
Meek,  (1980) ; Manson,  Meek  and  Gregory,  (1981)].  Saturation  of  gravity  wave  amplitudes 
occur  primarily  above  the  mesospheric  region  where  kinetic  viscosity  rises  above 
turbulence  viscosity  until  eventually  the  kinetic  viscosity  damps  out  the  gravity 
wave  mode  at  very  high  altitudes,  (above  500  km).  Since  some  important  infrared 
emissions  such  as  the  OH  emission  peaks  in  the  mesospheric  region,  the  instability 
can  produce  large  spatial  and  temporal  variations  in  infrared  radiance.  For 
part  (4),  we  investigate,  to  lowest  order,  the  production  of  a  secondary  short- 
period  oscillation  from  a  primary  long-period  gravity  wave.  The  instability 
can  produce  turbulence  and  secondary  oscillations  closely  associated  with  the 
Brunt  oscillation  with  frequencies  and  wave  lengths  which  may  be  analyzed  by 
the  dispersion  relation  derived  in  part  (3). 


Part  (1) 


(a)  The  Gravity  Wave  Model 


In  this  part  we  consider  the  effect  of  gravity  waves  on  the  number  density 
profiles  of  minor  atmospheric  constituents.  We  assume  that  the  gravity  waves 
are  carried  by  the  major  atmospheric  constituents  and  that  the  ambient  atmosphere 
(including  both  major  and  minor  components)  is  uniform  along  the  horizontal 
directions  (horizontal  stratification).  Using  the  theory  developed  by  Tuan  (1979) 
and  Hatfield  et  al  (1981),  we  can  then  show  that 
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concentration  of  the  ambient  minor  atmospheric  constituent 
perturbation  amplitude  of  the  minor  atmospheric  constituent 
horizontal  wave  vector  of  the  gravity  wave 
horizontal  particle  velocity  of  the  gravity  wave 
vertical  particle  velocity  of  the  gravity  wave, 
frequency  of  the  gravity  wave. 


The  assumption  of  horizontal  stratification  for  all  perturbed  physical  quantities 
(i.e.  AN,  Ap,  etc.)  means  that  they  may  take  on  the  generic  form: 
i(wt  -  k  x) 

f (z)e  X  (2) 


where  x  is  a  horizontal  co-ordinate.  Furthermore,  the  gravity  wave  must  satisfy 

Aw  =  o  at  the  ground  level,  (the  rigid  surface  boundary  conditions).  This  means 
8 

that  we  have  a  stationary  wave  along  the  vertical  direction;  and  that  all  the 
perturbed  quantities  must  also  be  stationary  along  the  vertical.  Thus,  f (z)  is 
either  pure  real  or  pure  imaginary.  The  imaginary  number  "i"  in  equation  (1) 
comes  from  the  complex  notation  we  use  and  provides  for  the  phase  difference 


-  5  - 


between  the  first  and  the  second  terms. 


For  the  gravity  waves  we  use  the  model  developed  by  Tuan  (1976)  and  Yu  et 
al  (1980).  The  model  expresses  the  horizontal  and  vertical  velocity  fields  in 
terms  of  the  pressure  perturbation,  Ap,  as  follows: 
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(0^  =  (^(z)  =  Brunt  frequency 


p  =  ambient  atmospheric  density. 


c'  =  derivative  of  the  speed  of  sound  profile. 


The  rigid  surface  boundary  condition  requires  that  Aw 
corresponding  condition  on  \p  is  given  by 
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In  this  way,  the  horizontal  and  vertical  velocity  fields  (Au  and  Aw  )  can 

g  g 

be  both  expressed  in  terms  of  the  pressure  perturbation  \p.  It  can  be  readily 
shown  (Tuan  (1976),  Yu  et  al  (1980)  that  satisfied  the  following  second  order 


self-adjoint  hydrodynamic  equation  : 
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+  g(  — j  )  =  Brunt  frequency. 


Equation  (6)  is  identical  to  the  usual  coupled  hydrodynamic  equations  in 
the  absence  of  dissipation.  It  is  possible  through  elimination  of  the  other 
variables  such  as  the  density  Ap,  the  horizontal  velocity  Au  and  the  vertical 
velocity  Aw  to  obtain  equation  (6)  for  ip  which  is  proportional  to  the  pressure 

O 

perturbation  Ap. 


Although  both  Q  and  V  in  equation  (6)  depend  on  the  frequency,  for  any  gravity 


wave  with  period  greater  than  30  minutes  Q  and  V  become  essentially  independent 


of  frequency.  Since  most  gravity  waves  found  in  nature  tend  to  have  periods 
of  the  order  of  one  to  two  hours  (Battaner  &  Molina  (1980)),  we  may  neglect  the 


variations  with  frequency  of  Q  and  V. 

To  calculate  Q  and  V  Tuan  (1976)  and  Yu  et  al  (1980)  have  used  the  1972 
COSPAR  model  for  the  temperature  profile,  the  mean  molecular  mass,  the  ratio 

t 

of  specific  heats.  The  variation  in  the  acceleration  due  to  gravity  is  of  course 
very  easily  taken  care  of.  For  the  present  report  we  have  used  5  different  exospheric 
temperatures.  They  are  600°k,  1000°k,  1500°k,  1800°k  and  2200°k.  These  temperatures 
taken  from  COSPAR  (1972)  are  the  closest  approximations  we  can  find  to  cases 
A  (626. 5°k),  C(1081 . 2°k) ,  F(1555.7°k),  G(1829.4°k)  and  H(2176.8°k)  in  the  AFGL 
Report  (AFGL-TR-81-0207) . 

We  have  used  gravity  wave  models  with  different  horizontal  phase  velocities 

corresponding  to  different  vertical  wave  lengths.  In  general,  it  is  the  vertical 

wave  lengths  that  produce  the  spatial  variations  in  infrared  radiance,  but  the 

vertical  wave  length  varies  with  height  in  our  inhomogeneous  atmosphere.  So 

the  horizontal  phase  velocity,  V  ,  =  ,  whose  inverse  is  just  ft,  the  square 

p  x 

root  of  the  eigenvalue  in  equation  (6)  and  can  therefore  be  used  to  designate 
the  type  of  gravity  wave  we  use.  In  general,  a  smaller  vertical  scale  structure 
implies  a  smaller  vertical  wave  length  which  in  turn  implies  a  smaller  horizontal 
phase  velocity.  From  the  paper  by  Tuan  (1976)  and  Yu  et  al  (1980)  we  can  show 
that  there  are  two  types  of  gravity  waves.  One  type  has  a  large  horizontal  phase 
velocity  (greater  than  400  m  sec  *)  with  its  energy  primarily  concentrated  in 
the  F-region  (above  110  km).  Such  waves  do  not  significantly  effect  the  lower 
altitudes  where  some  of  the  more  important  infrared  emission  profiles  peak. 

The  second  type  has  a  low  horizontal  phase  velocity  (less  than  400  m  sec  *)  with 


much  of  its  energy  concentrated  in  the  lower  atmosphere.  It  is  such  waves  which 
have  a  strong  effect  on  infrared  emission  that  we  use  as  gravity  wave  models. 


The  particular  choices  of  gravity  wave  horizontal  phase  velocities  are  based 
on  the  frequency  of  occurrence  of  typical  gravity-wave  velocities.  Battaner 
and  Molina  (1980)  have  found  that  the  horizontal  phase  velocities  typically  range 
from  60  to  140  m  sec  Since  the  more  pronounced  vertical  scale  structures 
are  produced  by  gravity  waves  with  low  horizontal  phase  velocities,  we  choose 
two  horizontal  phase  velocities  for  our  gravity-wave  models  in  the  lower  end 
of  the  typical  range,  60  and  78  m  sec  * .  We  also  choose  a  horizontal  phase  velocity 
of  129  m  sec  *  for  case  F.  Our  gravity  wave  model  is  valid  up  to  120  km.  Above 
120  km,  dissipation  due  to  molecular  viscosity  and  heat  conduction  gradually 
become  significant  and  the  models  cannot  be  used  without  modification. 

(b)  Spatial  variations  of  minor  atmospheric  constituents 

In  this  section  we  consider  the  effects  of  gravity  waves  on  the  concentration 

profiles  of  the  following  atmospheric  constituents  associated  with  infrared 

emission:  (1)  Day  [0] ;  (2)  Night  [0];  (3)  Day  [NO];  (4)  Night  |N0];  (SHCC^l; 

(6)  [H^O] .  For  each  ambient  profile  we  computed  the  perturbation  amplitude  AN 

from  equation  (1)  in  which  we  substitute  each  of  the  above  six  profiles  for  NQ  in 

turn.  Au  and  Aw  ,  the  gravity  wave  velocities  are  computed  from  equation  (3), 
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(4)  and  (6)  for  the  exospheric  temperatures  mentioned  in  part  (a)  which  correspond 
to  the  5  cases  A,  C,  F,  G  and  H  in  the  AFGL  Report.  Table  I  provides  the 
key  to  the  graphs  for  the  ambient  and  gravity-wave  perturbed  profiles. 

Fig.  9  to  38  (inclusive)  show  the  calculations  for  Day  [0],  Night  [0],  Day 
[NO],  Night  [NO] ,  [CO2  ]  and  [1^0]  profiles  under  the  influence  of  gravity  waves 
with  horizontal  velocities  of  60  to  78  m  sec  For  such  horizontal  phase  velocities, 
the  vertical  wave  length  is  small  compared  with  the  structure  scale  sizes  of 
the  unperturbed  atmospheric  constituent.  The  spatial  variations  in  the  profile 
imposed  by  the  gravity  wave  show  up  rather  dramatically.  This  would  not  be  true 


TABLE  I 
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for  a  gravity  wave  with  large  horizontal  phase  velocity  and  hence  greater  vertical 
wave  length.  If  the  vertical  wave  length  becomes  large  compared  with  the  scale 
size  of  the  background  atmospheric  structure,  no  spatial  variations  can  be  detected 
as  we  will  show  in  the  next  section. 

In  Figure  4,  5,  6,  7  and  8  we  examine  what  happens  when  an  intermediate 
horizontal  phase  velocity  (Vp^x  =  129  m  sec  *)  with  an  effective  larger  vertical 
wave  length  is  used.  The  results  show  some  vertical  structures.  The  structure 
produced  in  the  CO2  profile  (Fig.  2)  occurs  at  about  110  km  and  bears  some 
resemblance  to  the  SPIRE  data  (14.05  -  16.0  um)  •  Owing  to  the  natural  growth 
of  wave  amplitudes  with  height,  there  should  be  more  prominent  structures  at 
higher  altitudes.  However,  due  to  the  rapid  increase  in  molecular  viscosity 
above  110  km,  there  should  be  considerable  damping  of  wave  amplitudes  above  120 
km  leaving  the  most  prominent  oscillation  at  almost  110  km.  But  even  with  all 
the  dissipation  included,  it  is  difficult  to  see  how  all  traces  of  oscillation 


can  immediately  cease  above  110  km  as  the  SPIRE  data  seems  to  show.  In  fact, 
the  damping  of  a  gravity  wave  does  not  become  total  until  above  500  km.  Thus, 
in  the  absence  of  further  experimental  evidence,  we  feel  we  may  continue  to  support 
our  previous  conclusion  given  in  the  "Assessment  on  CO2  (V2)  Radiance  Profile," 
by  Chow  and  Tuan  (1980). 

(c)  Temporal  variations 

In  order  to  investigate  the  temporal  variability  of  infrared  radiance  produced 
by  gravity  waves,  we  have  considered  (1)  the  case  in  which  the  vertical  gravity 
wave  length  is  less  than  the  background  atmospheric  structure  and  simultaneous 
spatial  and  temporal  variations  may  be  observed;  (2)  the  cases  in  which  the  vertical 
wave  length  is  long  compared  with  background  structure  and  as  a  result  spatial 
variations  can  no  longer  be  easily  detected  but  temporal  variations  continue 
to  be  very  readily  observable.  In  (2)  we  consider  the  special  case  of  the  guided 
Lamb  mode  with  a  horizontal  phase  velocity  of  over  300  m  sec  *  and  a  rather  long 
vertical  wave  length.  The  choice  of  periods  for  the  gravity  waves  is  dictated 
by  typical  values.  Most  frequently  observed  gravity  waves  have  a  period  of  between 
2  and  3.5  hours  (Battaner  et  al,  (1977)).  We  have  chosen  a  gravity  wave  with 
a  period  of  2  hours. 

The  results  for  (1)  mentioned  above  are  shown  in  Fig.  la,  lb,  lc  and  Id  where 
profiles  for  Day  [0]  are  given  for  T  ■  0,  0.25,  0.5,  0.75  and  1  hour  respectively 
for  a  gravity  wave  with  a  2  hour  period  and  a  horizontal  phase  velocity  of  78 
m  sec  * .  Since  the  vertical  wave  length  for  such  a  wave  is  small  compared  with 
the  scale  size  of  the  background  [0]  structure,  we  see  considerable  simultaneous 
spatial  and  temporal  variations.  In  Fig.  2  we  see  what  happens  to  the  same  profile 
when  the  horizontal  phase  velocity  is  close  to  a  guided  mode,  in  particular, 
the  Lamb  mode.  For  such  a  mode  the  vertical  wave  length  is  long  compared  to 


the  scale  size  of  the  ambient  profile  structure  and  we  see  that  the  profile  shape 
is  unaffected  by  the  wave.  Instead,  the  whole  profile  oscillates  up  and  down 
at  T  -  0  and  1  hour  for  a  gravity  wave  with  a  2-hour  period.  The  only  infrared 
variance  associated  with  atomic  oxygen  that  may  be  observed  will  be  through  the 
ground-based  photometers  which  would  show  a  periodic  oscillation  in  the  total 
zenith  intensity.  This  would  be  the  results  for  (2)  mentioned  above  where  spatial 
variations  is  all  but  unobservable  and  one  can  only  observe  the  temperal  changes. 
Fig.  3  shows  the  same  type  of  phenomena  for  Day  [NO]  profile 


Part  (2) 

The  infrared  volume  emission  rate  profiles  are  calculated  from  the  concentration 
profiles  mentioned  in  part  ( 1  )  by  assuming  that  all  such  emissions  come  from 

i 

I 

,  i 

spontaneous  emission  of  vibrationally  excited  states  of  [0],  [NO],  [C0?]  and 
TH^O]  (See  AFGL  Report).  The  problem  therefore  reduces  to  the  determination 
of  the  concentration  of  vibrationally  excited  states.  For  some  infrared  emitters 
such  as  [OH]  the  principal  excitation  mechanism  is  through  its  chemical  formation 
(H  +  0^  OH*  +  0^)  which  leaves  it  in  an  excited  vibrational  state.  For  the 
minor  atmospheric  constituents  we  treat  here  in  this  report,  the  principal  excitation 
mechanism  is  through  absorption  of  infrared  radiation.  Such  absorption  can  come 
from  three  different  radiation  sources:  (1)  the  radiation  from  the  lower 
atmospheres  considered  as  a  plane  blackbody;  (2)  the  direct  solar  radiation; 

(3)  the  radiation  from  other  atoms  and  molecules  in  the  atmosphere.  Since  all 
three  radiation  sources  have  to  pass  through  a  certain  thickness  of  atmosphere 
before  they  reach  the  height  level  z  where  we  compute  the  concentration  of  excited 
states,  we  need  to  use  the  well-known  radiative  transfer  functions  to  take  care  4 

of  the  loss  of  intensity  of  the  radiations  as  they  traverse  from  their  source 
through  the  atmosphere  to  reach  the  point  at  z.  Such  radiative  transfer  functions 


are  not  necessary  for  OH  infrared  emission  since  its  principal  excitation  mechanism 
is  through  local  chemical  reactions  rather  than  through  radiation  sources  from 
far  away. 

Once  the  excitation  rate  is  computed,  the  concentration  of  each  particular 

excited  state  may  be  determined  through  assuming  chemical  equilibrium  (AFGL  Report) . 

Thus,  if  [X]  is  the  concentration  of  the  minor  species  which  emit  infrared  radiation 
* 

and  [X]  is  the  concentration  of  a  particular  vibrationally  excited  state  and 

* 

[Y]  and  [Z]  are  concentrations  of  any  two  gases  which  collide  with  [X]  and  [X] 
respectively,  we  may  write  the  equation  for  chemical  equilibrium 


ot[X]  [Y]  +  B[X]  -  6fX]*[Y]  +  AfX]*  (7) 

where  a  ■  rate  coefficient  for  collisional  excitations 

B  =  radiative  excitation  coefficient  and  is  calculated  from  the 
radiative  transfer  functions. 

8  ■  rate  coefficient  for  collisional  de-excitation 
A  <=  Einstein  coefficient  for  spontaneous  radiation  of  [X]  . 

Equation  (7)  essentially  assumes  that  production  rate  of  excited  states  is 
equal  to  the  loss  rate  and  is  sufficiently  accurate  at  altitudes  up  to  200  km. 

It  is  an  algebraic  equation  from  which  the  concentration  of  excited  states  [X] 
can  be  readily  solved. 

At  higher  altitudes,  however,  there  is  an  additional  loss  rate  through  vertical 

diffusion  which  is  neglected  in  the  calculations  in  AFGL  Report  (AFGL-TR-81- 

0207)  and  which  may  be  important.  To  take  this  into  account  we  have  to  write 

* 

down  the  following  mass  and  momentum  conservation  equation  for  [X]  : 


am* 


afXHY]  +  B[X]  -  8[X]*[Z1  -  AfXl* 


-  div  {[X]  vd  } 
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m[X]  g 


(9) 


[X]*  kT 


where  v  =  collision  frequency  with  other  molecules 


p  =  partial  pressure  of  [X] 


=  vertical  drift  velocity 


X 

In  equation  (8),  we  may  assume  quasi-equilibrium  (i.e.  —  «  0).  This 

assumption  was  also  made  in  the  AFGL  Report.  Prom  equations  (9)  and  (10)  we  obtain 
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where  D  =  —  *  the  diffusion  coefficient  for  [X] 
mv 


H  =  —  =  scale  height  for  [X] 

mg 

In  deriving  (11)  we  have  assumed  that  the  kinetic  temperature  in  nearly  constant 
above  200  km  where  the  diffusion  term  becomes  important.  Substituting  (11) 
in  (8)  and  setting  3fX]  /3t  ■  0  ,  we  obtain 


g|Df  gf  +  [X]*  -  (A  +  8[z))[XT*  +  <x[X][Y]  +  B[X]  -  o 


Quite  obviously  equation  (12)  reduces  to  equation  (7),  the  equation  for 


chemical  equilibrium  when  we  neglect  the  diffusion  (i.e.  setting  D  -  0).  Unlike 

(7)  equation  (12)  is  a  second  order  self-adjoint  differential  equation  and  requires 

h 

more  work  to  obtain  [X]  ,  the  concentration  of  the  excited  states. 

The  physical  reason  why  the  diffusion  term  is  important  is  that  the  diffusion 

coefficient  is  inversely  proportional  to  the  collision  frequency.  At  higher 

altitudes,  there  will  be  a  rapid  decrease  in  the  collision  frequency,  thus  producing 

a  large  Increase  in  D(z).  Such  an  increase  in  vertical  diffusion  will  result 

* 

in  an  increased  in  [X]  at  higher  altitudes. 

Fig.  39  to  44  (inclusive)  show  the  effect  of  gravity  waves  with  a  horizontal 
phase  velocity  of  78  m  sec  *  on  both  the  day  time  and  the  night  time  infrared 
emission  rate  profiles  of  [NO] .  We  have  shown  only  the  1  0  vibrational  transition 

with  a  band  center  at  5.325  ym.  It  is  interesting  to  compare  these  results  with 
Fig.  11,  12,  21,  22,  31  and  32  which  show  the  effect  of  gravity  waves  on  the 
concentration  profiles.  It  is  obvious  that  the  variations  are  greater  in  the 
concentration  profiles  than  they  are  in  the  emission  profiles.  The  physical 
reasons  for  this  are  easily  understood.  For  [0],  [NO],  [CO^]  and  [H^O]  the  excited 
states  are  produced  principally  by  the  absorption  of  radiation  from  three  different 
sources;  the  infrared  radiation  form  the  lower  atmosphere  considered  as  a  plane 
blackbody;  the  direct  solar  radiation;  the  radiation  from  the  rest  of  the  atmosphere. 
The  latter  effect  has  to  be  integrated  over  the  entire  height  range  of  the  atmosphere 
and  has  a  significant  smoothing  effect. 

There  are  other  types  of  infrared  emissions  such  as  OH  emissions,  for  instance, 
in  which  the  principal  excitation  mechanism  is  through  local  chemical  reaction 
(i.e.  H  +  0^  OH  +  O2)  and  there  is  no  "smoothing"  of  the  emissions  profile 
from  an  integration  over  the  height  range  of  the  atmosphere.  The  production 
of  excited  OH  states  are  then  completely  local  and  the  emission  profiles  can 
exhibit  far  more  pronounced  spatial  structure.  This  can  be  seen  from  the  computed 


profiles  of  Hatfield  et  al  (1981)  in  which  the  magnitude  of  the  variations  in 
the  spatial  structures  of  the  emission  profiles  are  comparable  to  the  corresponding 
variations  in  the  concentration  profiles. 


Part  III 

To  analyze  the  more  sharply  pronounced  spatial  and  temporal  structure  with 
scale  sizes  of  the  order  of  a  few  tens  of  kilometers  and  a  few  minutes,  we  need 
to  develop  a  dispersion  relation  for  the  short  period  gravity  waves  with  frequencies 
close  to  the  Brunt  frequency.  Since  such  spatial  and  temporal  inhomogenieties 
tend  to  occur  in  the  mesospheric  region  where  variations  in  the  temperature  and 
mean  molecular  mass  of  the  ambient  atmosphere  as  well  as  the  possibility  of  heavy 
wave  reflection  from  the  sharp  temperature  rise  in  the  base  of  the  thermosphere 
require  drastic  modifications  to  the  Hines'  dispersion  relation  (Hines,  (1960); 

Eckart,  (I960)).  Actually,  generalizations  of  this  dispersion  relation  to  an 
inhomogeneous  but  horizontally  stratified  atmosphere  including  dissipation  has 
been  made  in  thw  WKB  approximation,  (Pitteway  &  Hines  (1963),  Midgley  and  Liemohn 
(1966),  Yeh  and  Liu  (1974),  Francis  (1975)). 

It  is  the  purpose  of  the  present  paper  to  analyze  a  generalization  of  Hines' 
dispersion  relation  to  an  inviscid  but  inhomogeneous  and  horizontally  stratified 
atmosphere  which  remains  valid  not  only  for  free  modes,  but  also  for  guided  or 
partially  guided  modes.  In  the  derivation  we  use  methods  and  techniques  based 
on  a  "potential"  treatment  of  gravity  waves  (Tuan  (1976),  Yu  et  al  (1980),  Tuan 
and  Tadic  (1982)).  The  generalized  dispersion  relation  would  retain  a  relatively 
simple  analytic  form  and  can  hence  be  easily  used  to  estimate  wave  lengths,  frequencies 
or  horizontal  phase  velocities  as  the  case  may  be  as  well  as  to  provide  simple 
qualitative  estimates  of  the  limits  of  validity  of  the  Hines'  dispersion  relation. 

To  accomplish  this,  we  begin  by  examining  how  a  vertical  wave  vector  and  hence 
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a  "local"  vertical  wave  length  can  be  meaningfully  defined  In  an  Inhomogeneous 
atmosphere.  Two  specific  cases  will  be  studied.  One  case  will  deal  with  gravity 
waves  with  frequencies  close  to  the  Brunt  and  acoustic  cut-off  frequency  while 
the  other  case  deals  with  long  period  gravity  waves.  For  the  former  case  we 
will  show  that  in  altitude  rangesvhere  tu^Cz)  "*■  w,  the  "local"  vertical  wave  vector 
will  always  be  very  large  and  pure  imaginary.  For  the  latter  case  we  will  examine 
the  variation  of  horizontal  phase  velocity  with  height  for  zero  vertical  wave 
vector.  From  such  a  curve,  one  can  identify  a  region  of  horizontal  phase  velocities 
in  which  the  Hines'  relation  is  valid.  Such  a  curve  will  also  show  for  a  fixed 
horizontal  phase  velocity  the  height  ranges  in  which  the  Hines'  dispersion  relation 
is  valid. 

In  the  appendix  we  will  work  out  the  lower  frequency  limit  as  well  as  the 
height  range  in  which  the  generalized  dispersion  relation  can  be  be  used. 


In  an  inviscid  inhomogeneous  and  horizontally  stratified  atmosphere,  the 
linearized  hydrodynamic  equations  for  monochromatic  wave  solutions  are  given  by 


iuAp  “ 

ik  p  Au 

X  0 

-  Ji  (0o4v) 

(13) 

iwp  Au 

0 

=  ik  Ap 

X 

iupoAw 

=  -gAp  - 

3tfp) 

9z 

iwAp  = 

iu>  . 

— 2  Ap  + 

✓—V 

N 

V- ' 

CM 

3°  “ 
o 

Q. 

where  Ap,  Au,  Aw  and  Ap  are  respectively  the  density,  horizontal  velocity,  vertical 
velocity  and  pressure  perturbations,  w.(z)  is  the  Brunt  frequency,  w  is  the 


gravity  wave  frequency,  is  the  horizontal  wave  vector,  and  c  is  the  speed 

of  sound  profile.  Tuan  and  Tadic  (1982)  have  shown  that  equation  (13)  can  be 
transformed  into  the  following  Sturm-Liouville  equation  for  the  pressure  perturbation 
without  approximations. 
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a  =  some  upper  boundary  above  which  the  atmosphere  may  be  assumed  to  be  uniform. 
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In  equation  (14),  V (a)  ■  0  and  W(a)  *  1.  Thus,  above  the  upper  boundary, 

$  has  only  plane  wave  solutions.  Below  the  upper  boundary  V(z)  and  W(z)  are 
well-behaved  functions  of  z  as  well  as  the  parameter  0),  (Their  explicit  forms 
are  given  by  eqution  (10)  in  the  paper  by  Tuan  and  Tadic  (1982)). 

Equation  (13)  may  be  put  into  the  form: 


[  ~K  +  K  2  (z)  ]  *  =  0 

dz2  _z 

In  this  form  0  may  be  considered  as  a  solution  of  a  wave  equation  with  local 
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wave  vector  K  and  local  wave  length  —  • 

z  K 

— z 

To  see  this  from  another  point  of  view  we  can  divide  the  atmosphere  up  into 
thin  uniform  horizontal  slabs  such  that  V  and  W  are  constant  within  each  slab. 


This  means  of  course  that  K  (z)  is  also  constant  within  each  slab.  Thus,  the 

— z 

solutions  in  the  i  th  slab  with  thickness  Az  are  given  by: 


m 

1$ 

m 


* 


£5 


<r.  •: . 


tf3 


>.\  •  ■ 


S3 


These  solutions  in  (17)  are  obviously  plane  wave  solutions  with  vertical  wave 
2n 


length  X 
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We  can  construct  the  entire  solution  for  all  z  by  imposing 


boundary  conditions  at  the  ground  and  the  upper  boundary  as  well  as  suitable 
interfacial  boundary  conditions  between  slabs.  In  this  way  the  local  wave  length 
X^_  may  be  dfefined  by  shrinking  the  slab  to  very  small  thickness  so  that 
X  =  X  .. 

Z  A  n  zi 

Az^  -*■  0 

This  approach  is  quite  different  from  the  usual  multi-layer  method  such  as 
the  one  used  by  Volland  (1969)  in  that  the  slabs  are  not  necessarily  uniform 
in  temperature  or  mean  molecular  mass  but  are  uniform  in  the  quantities  V(z) 
and  W(z)  which  have  been  previously  identified  as  "potentials"  in  analogy  with  the 
Schrodlnger  equation  (Yu  et  al  (1980),  Tuan  and  Tadic  (1982)).  Thus,  instead 
of  isothermal  slabs,  we  have  essentially  "isopotential"  slabs.  The  difference 
between  the  isopotential  and  the  isothermal  slab  is  that  all  height  derivatives 
of  atmospheric  parameters  (such  as  temperature,  background  density  etc)  are 
automatically  included  in  the  isopotential  slabs. 

Here  we  wish  to  emphasize  that  the  local  wave  length  we  have  used  is  exact 
and  not  derived  through  a  WKB  type  of  approximation. 

By  substituting  the  expressions  for  V(z)  and  W(z)  (taken  from  equation  (10), 
Tuan  and  Tadic  (1982))  we  can  show  that 
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(  )2  +  g(  — ^  )  =  acoustic  cut-off  frequency; 


and  all  primes  on  the  speed  of  sound  profile  c(z)  denote  height  derivatives. 

It  is  obvious  that  for  an  uniform  isothermal  atmosphere  =  ^1  ~  ^2  “  ® 

and  equation  (18)  immediately  reduces  to  the  Hines'  dispersion  relation  given 
by  the  bracket  term  on  the  right-hand  side  of  (18).  It  Is  interesting  to  note 
that  since  £q ,  and  ^  depend  only  on  the  background  atmosphere  and  do  not 
depend  on  the  gravity-wave  parameters,  at  low  frequencies  (u)  «  ui^)  ,  all  three 
correction  terms  to  the  Hines'  relation  become  independent  of  the  gravity-wave 
frequency.  Thus,  if  the  period  is  long,  the  Hines'  relation  is  true  only  for 
small  horizontal  phase  velocity  for  which  the  first  term  on  the  right-hand-side 
of  equation  (18)  dominates. 

For  short-period  gravity  waves  the  most  interesting  term  in  equation  (18)  is 
the  last  term  on  the  right-hand  side.  Quite  obviously,  for  any  horizontal  phase 
velocity,  as  -*  w  the  last  term  on  the  right-hand  side  of  (18)  dominates.  Sin< 
£0(z)  >  0  for  all  z  and  the  denominator  is  always  positive,  this  term  is  always 


negative.  This  also  means  that  at  those  height  levels  where  w^(z)  ~  w,  K is 

always  pure  imaginary  and  very  large.  Hence,  the  wave  is  always  purely  evanescent 

and  the  amplitude  drops  to  zero  very  quickly.  To  see  this  we  consider  the  "potential 

(Tuan  1976,  Yu  et  al  (1980))  for  a  gravity  wave  with  a  7.5  minute  period,  (Fig. 

45).  Quite  obviously  at  approximately  160  to  170  km  there  is  a  sharp  infinite 

2  2 

potential  barrier  produced  by  the  vanishing  of  the  denominator  (z)  -  u  .  Within 

this  barrier  is  infinite  and  pure  imaginary  and  there  is  no  wave  propagation 

2  2 

within  this  region.  Above  this  region  <  to  and  only  acoustic  waves  can  propagate 

The  infinite  potential  barrier  therefore  separates  the  gravity-wave  mode  from 

the  acoustic-wave  mode.  This  can  be  seen  in  Fig.  (46)  in  which  a  7.5  minute 

period  wave  would  cut  the  Brunt  period  T^  at  163  km  where  the  "potential"  rises 

to  infinity.  At  167  km  it  intersects  with  the  acoustic  cut-off  period  so  that 

above  167  km  acoustic  wave  propagation  is  permissible.  Since  we  have  identified 

K  as  the  vertical  wave  vector,  we  shall  from  now  on  refer  to  both  K  and  k 
— z  — z  z 

as  simply 


To  facilitate  the  use  of  the  dispersion  relation  (18),  we  have  provided  a 

2  2  2 

table  (Table  II)  for  the  function  ,  c  ,  £0»  and  ^  at  5  km  intervals 

from  the  ground  level  up  to  120  km.  Since  £Q*  ^  and  depend  on  derivatives 
of  atmospheric  parameters,  the  results  are  exceedingly  sensitive  to  the  model 
atmosphere.  We  have  used  the  Winter  mid-latitude  COSPAR  1972  model.  Great  care 
must  be  taken  in  smoothing  a  curve  for  further  numerical  differentiation,  since 
this  often  results  in  shifting  the  phase  of  the  differentiated  curve.  We  have 
adopted  the  inverse  procedure  of  smoothing  the  differentiated  curve  and  integrating 
to  verify  that  we  obtain  the  same  original  curve. 

We  compare  equation  (18)  with  the  Hines’  dispersion  relation  for  two  cases. 


One  for  w 


~  “b 


In  Fig.  46  there  are  three  vertical 


and  the  other  for  to  << 


lines  corresponding  to  gravity-wave  periods  T  =  3.75,  4.6  and  6.1 


minutes.  These  lines  either  come  close  to  or  intersect  with  the  Brunt  period 


at  different  altitudes.  At  such  altitudes  w  ~  ax  and  two  out  of  the  three  corrections 


to  the  Hines'  dispersion  relation  become  important.  Table  III  shows  the  values 


of  as  calculated  by  the  Hines'  relation  and  equations  (4)  at  the  above  altitudes. 


The  values  of  k  are  so  chosen  that  k  is  either  immediately  above  or  below 

X  2 


zero,  (i.e.  close  to  the  turning  point  of  the  WKB  approximation), 


In  general,  one  would  expect  the  Hines'  relation  as  well  as  the  WKB  approximation 


to  become  very  inaccurate  close  to  such  turning  points.  As  can  be  seen  from 


Table  III  there  is  an  enormous  difference  between  kz  obtained  from  (16)  and 


the  Hines'  relation. 


For  very  much  larger  gravity-wave  periods  (such  as  120  minutes)  we  may  neglect 
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w  in  comparison  with  or  and  the  corrections  to  the  Hines  '  relation  can 


continue  to  be  important  for  values  of  k^  close  to  zero  but  are  in  general  much 


smaller.  Table  III  again  shows  the  significant  albeit  much  smaller  difference 


between  the  values  of  k  obtained  from  (18)  and  from  the  Hines'  relation.  For 

z 


all  long  period  gravity  waves  the  correction  to  the  Hines'  relation  does  not 


arise  as  a  result  of  the  vanishing  of  the  denominator  (w^  -  u>  ) ,  but  as  a  result 


of  large  enough  horizontal  phase  velocities  for  which  the  Hines'  relation  in 


the  bracket  term  of  (18)  is  all  but  cancelled  by  the  three  correction  terms  which 


no  longer  depend  on  gravity-wave  parameters. 
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Fig.  47  shows  a  plot  of  V  (kz  =  0)  against  altitude  for  large  period 


2  2 

gravity  wave.  From  (18)  and  the  assumption  w  «  o>^  ,  we  obtain: 
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where  f (z)  = 


a)  2 3 4(z) 
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^(z)  C2(z) 
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“b  (V 


For  any  k  we  can  write 
z 


phx,o 


To  take  a  specific  example,  we  shall  assume  that  the  criterion  just  mentioned 
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is  satisfied  if  V  V  ^  o  *  ^len  we  just  simply  draw  the  curve 

2  12 

V  ,  =  -T7T  V  ,  and  we  know  that  in  the  shaded  region  below  this  curve  the 

phx  10  phx.o 

first  term  on  the  right-hand  side  of  (8)  which  is  the  important  term  in  the  Hines' 
relation  is  at  least  10  times  greater  than  the  "correction"  to  the  Hines'  relation. 
Hence,  we  may  well  expect  that  within  this  shaded  region,  the  Hines  relation 
is  a  good  approximation. 

2  1 

On  the  other  hand,  if  we  decide  that  for  V  ,  to  be  in  the  range  _  V  , 

phx  °  2  phx,o 

2  1  2 

>  Vp^x-  >  Jo  Vphx  o  ’  ^le  ^nes*  relation  is  good  enough,  we  can  draw  in  the 


curve  V 


1  2 

~2  Vphx  *  0  anc*  assume  that  within  this  region  the  Hines'  relation 


2  1  1 

is  say  fair  to  good.  For  V  .  >  V  . 

3  6  phx  2  phx,o 

not  particularly  good. 


The  Hines'  relation  is  obviously 


As  an  example  of  how  one  may  make  use  of  Fig.  47  we  consider  a  long-period 


gravity  wave  with  a  horizontal  phase  velocity  of  200  m  sec  .  We  draw  a  horizontal 


line  in  Fig.  47  for  V 


2  -2 

0.04  km  sec  .  We  see  that  between  30  km  and  77.5  km 


as  well  as  above  109  km  the  horizontal  phase  velocity  falls  within  the  region 

4  V  ,  2  >vt^>T^V,2  within  which  the  Hines'  relation  is  fair  to  good. 

2  phx,o  phx  10  phx 

It  should  not  be  used  at  the  other  height  ranges.  Here  we  should  emphasize  that 


the  numbers  1/2  and  1/10  are  chosen  purely  arbitrarily  to  illustrate  our  point. 
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Clearly  the  curve  V 


phx 


2  2 
Vphx  Q  separates  the  region  k^  >0  from  the 


region  k^  <0.  Thus,  above  this  curve  the  gravity  wave  becomes  evanescent. 

Our  dispersion  relation  given  by  (6)  remains  valid  for  all  these  regions.  We 
can  hence  use  Fig.  47  to  determine  the  regions  in  the  atmosphere  in  which  a  long- 
period  gravity  wave  with  a  given  horizontal  phase  velocity  may  be  evanescent. 

As  an  example,  we  consider  a  gravity  wave  with  a  horizontal  phase  velocity  of 

-1  2  2-2 
say  346  m  sec  .  This  means  V  ^  =  0.12  km  sec  .  The  horizontal  line  for 

this  phase  velocity  in  Fig.  47  shows  that  the  wave  is  evanescent  everywhere  except 

in  a  narrow  region  between  44  and  64  km  .  We  see  also  that  the  wave  can  propagate 

above  116  km. 


Part  (4) 

Whilst  a  large  variety  of  instabilities  may  occur  in  the  atmosphere,  it  has 
long  been  known  that  an  important  source  of  atmospheric  instability  and  turbulence 
in  the  mesospheric  region  may  be  caused  by  the  growth  of  gravity-wave  velocity 
amplitudes  with  height  until  Richardson's  number  falls  below  some  critical  value 
[Hodges,  1967;  Hines,  1974;  Roper,  1977;  Zimmerman  and  Murphy,  1977;  Manson  and 
Meek,  1980;  Manson  et  al  (1981)].  It  is  also  well  known  that  important  optical 
emissions,  such  as  the  5577A  01  and  the  OH  infrared  emissions,  also  peak  in  the 
same  region  (i.e.  between  80  km  and  100  km). 

Ground-based  photometric  measurements  have  shown  that  in  the  presence  of 
a  gravity  wave  the  5577A  01  emission  often  exhibits  a  short  period  oscillation 
(of  the  order  of  minutes)  superimposed  on  a  long  period  gravity-wave  oscillation  TOkuda 
1962;  Silverman,  1962;  Tuan  et  al. ,  1979].  OH  airglow  photographs  by  Peterson 
and  Kieffaber  from  1972  onward  show  an  extremely  wide  variety  of  moving  and 
stationary  airglow  structures,  some  of  which  vary  in  brightness  [Peterson  and 
Kieffaber,  1973;  Peterson,  1979].  Parallax  measurements  from  the  ground  and 
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from  aircraft  place  these  features  at  heights  of  80  to  100  km  [Peterson  and 
Kieffaber,  1973;  1975].  Of  particular  interest  are  airglow  features  of  extremely 
short  wavelength  and  short  time  duration  which  we  have  called  "ripples".  These 
"ripples"  (1)  have  a  very  uniform  wavelength  of  the  order  of  ten  kilometers; 

(2)  have  a  constant  horizontal  velocity  which  can  have  values  from  0  to  100  m/s; 

(3)  have  periods  of  the  order  of  a  few  minutes;  (4)  appear  suddenly  and  fade  with 
a  lifetime  of  about  45  minutes;  and  (5)  sometimes  appear  to  be  at  a  different 
height  from,  and  exhibit  relative  motion  with  respect  to,  the  large-scale  background 
airglow  structures.  In  addition  we  have  found  that  the  time  of  occurrence  of 
ripples  is  very  well  correlated  with  lower  lunar  transit  (LLT) .  Fig.  49  shows 

this  correlation. 

During  her  thesis  work,  Kieffaber  [1973]  made  many  all-sky  maps  at  1.65  u 
and  2.15  u  while  Peterson  photographed  the  sky  in  the  near  infrared.  Fig.  48 
shows  an  all-sky  map  for  the  night  of  2  December  1972.  Darker  areas  represent 
higher  intensity  of  OH  emission.  On  this  night  ripples  were  photographed  in 
the  brightest  region  of  OH  emission,  in  the  southwest  part  of  the  sky.  Lunar 
transit  occurred  at  2105  MST,  near  the  middle  of  the  mapping  period,  and  the 
ripples  were  photographed  at  2130  MST. 

In  general,  gravity-wave  oscillations  observed  through  optical  emissions 
are  divided  into  two  types.  One  type  shows  a  relatively  smooth  oscillation  with 
a  period  of  the  order  of  hours  ,  [Dachs,  1968].  The  second  type  shows  a  long 
period  oscillation  with  a  short  period  oscillation  of  the  order  of  minutes  superimposed. 
Figures  50  and  51  are  examples  of  these  second  types  of  waves.  The  small  period 
oscillations  usually  occur  only  in  the  presence  of  large  amplitude  gravity  waves. 

The  superimposed  short  period  oscillation  may  appear  gradually  as  a  self-excitation, 
or  it  may  appear  more  violently. 

Battaner  &  Molina  (1980)  have  used  the  first  type  of  waves  as  seen  in  the 
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5577  01  emission  for  an  estimate  of  the  turbulence  diffusion  coefficient.  The 
estimate  is  based  on  power-spectrum  analyses  of  the  observed  long-period  oscillation 
in  the  green-line  emission  coupled  with  the  theoretical  formula  for  the  turbulent 
diffusion  tensor  deduced  from  a  statistical  theory  of  turbulence  (Weinstock, 

1976).  The  results  appear  to  be  very  close  to  those  obtained  from  a  linearized 
monochromatic  theory  (Hines,  1960)  showing  that  most  gravity  waves  at  this  height 
level  (~90  km)  experience,  a  small  enough  frequency  spread  for  a  linearized  theory 
to  remain  valid. 

The  primary  objective  of  part  (3)  is  to  use  an  extension  of  the  «ory  of 
Tuan  et  al  (1979)  as  well  as  the  dispersion  relation  developed  in  part  (3)  to 
show  that  large  amplitude  primary  wave  oscillations  can  producte  secondary  oscillations 
with  periods  of  the  order  of  minutes  consistent  with  the  observed  temporal  variations 
in  the  5577A  green  line  as  well  as  the  spatial  variations  observed  by  Peterson. 

In  the  case  of  the  green  line,  the  small  period  oscillations  are  usually  observed 
only  when  the  primary  wave  has  rather  large  amplitudes  as  mentioned  in  the  introduction^ 
In  the  case  of  "ripples"  observed  by  Peterson,  whilst  a  fairly  significant  primary 
oscillation  with  a  period  of  about  2  hours  exists  on  the  night  of  2-4  Feburary 
1973,  (Fig.  54),  they  have  not  consistently  been  observed.  The  possibility  therefore 
remains  that  some  of  the  "ripples"  may  either  be  caused  by  other  mechanisms  or 
may  even  be  themselves  primary  gravity  waves.  In  this  part  of  the  report  we 
shall  be  only  concerned  with  gravity-wave  induced  instabilities. 

For  the  analysis  of  temporal  behaviour  of  the  small  period  oscillations, 

we  use  an  extension  of  the  theory  of  Tuan  et  al  (1979)  to  include  the  superadiabatic 
2  s  dO 

condition  (i.e.  =  "■  —  <  0  where  is  the  Brunt  frequency  and  0  is  the 

temperature  potential)  for  a  time  duration  short  compared  to  the  wave  period. 

The  latter  condition  can  occur  for  a  small  portion  of  a  gravity-wave  cycle  when 
its  amplitude  exceeds  some  critical  value  as  shown  in  Fig.  53.  The  amplitude 
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of  the  small  period  oscillation  following  a  short  duration  of  superadiabatic 
condition  can,  unlike  the  case  of  self  excitation  in  Fig.  52,  immediately  acquire 
relatively  large  values,  sometimes  comparable  to  the  amplitude  of  the  gravity 
wave  itself  as  in  Fig.  50.  The  mathematical  formulation  for  both  the  solutions 
shown  in  Fig.  52  and  Fig.  53  are  based  on  the  Semi-Lagrangian  formulation  of 
hydrodynamics.  In  this  formulation,  instead  of  looking  for  field  solutions, 
iKx\t)  where  \p  may  be  the  pressure,  the  density  or  the  velocity  fields,  we  consider 
the  trajectory  3?  (xq  ,t)  of  a  fluid  particle.  Since  we  may  in  general  neglect 
the  horizontal  motion  relative  to  the  vertical  owing  to  the  large  horizontal 
wave  length,  we  may  restrict  our  considerations  to  the  vertical  trajectory  z=z(zQ,t) 
By  substituting  such  a  solution  into  $  so  that  ^  =  lKz(t),t)  is  a  function  of 
only  the  time  variable  we  may  obtain  a  fairly  simple  ordinary  equation  for  z  =  z(t) 
if  we  neglect  higher  order  terms.  Thus,  by  substituting  =  ip(z(t),t)  into  the 
hydrodynamic  equations  and  expanding  the  solution  z(t)  about  z_(t)  =  X  sin  ut  +  z 

Vj  —O  O 

where  z  is  the  trajectory  of  a  fluid  particle  for  the  linearized  gravity  waves  and  i 

U 

is  the  primary  gravity  wave  frequency,  we  obtain  the  following  equation  for 
z  =  z(z  ,t)  to  lowest  order  in  (z  -z_). 

O  \j 


d  (z-zG)  d(z  -  zG) 

- __  +  n 


dt 


dt 


+  My  (t)(z  -  zG)  *  f(zG,t) 


(21) 


In  equation  (21) 

z G  *  sin  (t  -  tQ)  +  zQ  is  the  vertical  displacement  of  the  primary 
gravity  wave. 

=  a  coefficient  which  will  depend  on  the  turbulence  viscosity  in 
some  complicated  way  but  will  be  used  primarily  as  a  parameter 
for  this  report. 


n 


~  Wg  {[1  +  kzX  sin  [w(t  -  tQ)  -  kz(zG  -  zq)]} 


f(z  ,t)  =  forcing  functions  and  is  usually  small. 

Vjr 

We  see  immediately  that  the  superadiabatic  condition  occurs  if  k  X  >1 

z  o 

2 

when  turns  negative  for  a  portion  of  the  gravity-wave  cycle.  During  this 
period,  turbulence  is  generated,  causing  the  more  violent  small  period  oscillations 
to  follow,  Fig.  53. 

If  the  short-period  oscillations  should  give  rise  to  secondary  waves  the 
spatial  variations  of  such  waves  may  be  analyzed  by  the  dispersion  relation  discussed 
in  part  (3).  One  possible  mechanism  for  producing  such  waves  is  the  growth  of 
a  primary  gravity  wave  with  height  until  instability  occurs  at  some  regions  in 
the  mesosphere.  For  gravity  waves  guided  from  a  relatively  distant  source,  the 
wave  will  be  stationary  along  the  vertical  direction  and  instability  of  the  type 
just  discussed  can  occur  over  the  shaded  regions  shown  In  Fig.  54.  Such  instability 
can  occur  periodically  over  the  periods  of  the  primary  gravity  wave,  provided 
it  maintains  its  strength.  If  the  primary  gravity  wave  should  be  a  travelling 
wave,  then  the  region  of  instability  and  turbulence  would  propagate  along  with 
the  gravity  wave. 
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In  this  section  we  summarize  the  basic  physics  that  are  pertinent  to  spatial 
and  temporal  infrared  variance  produced  by  atmospheric  gravity  waves.  There 
are  in  general  two  types  of  gravity  waves.  The  type  which  produce  the  large- 
scale  T.I.D.'s  are  characterized  by  (see  Thome,  (1968); 

(i)  Energy  concentrated  at  high  altitudes  with  little  wave  amplitude 
below  110  km. 

(ii)  High  horizontal  phase  velocity  (in  excess  of  400  m  sec  *)• 

(iii)  Source  also  located  at  high  altitudes,  usually  from  polar 
magneto spheric  substorms. 

(iv)  Usually  travel  along  longitudes  and  exhibit  little  inclination 
to  go  past  the  equator. 

The  type  which  produce  the  mediun-scale  T.I.D.'s  are  characterized  by  (Thome, 
(1968),  Battaner  &  Molina  (1980),  Hines,  (I960)). 

(i)  Energy  mostly  concentrated  at  low  altitudes  but  can  spread  up  to  high 
altitudes 

(ii)  Low  horizontal  phase  velocities  (less  than  400  m  sec  *)• 

(iii)  Source  is  usually  at  low  altitudes  and  are  often  due  to 
meteorological  conditions 

(iv)  Do  not  travel  in  any  preferred  longitudinal  or  latitudinal 
direction. 

Gravity  waves  should  also  be  divided  into  two  categories: 

(i)  Those  with  small  enough  amplitudes  so  that  the  waves  may  be  accurately 
considered  as  purely  linear  waves.  These  gravity  waves  constitute  a  majority 

of  the  waves  which  are  far  from  their  sources. 

(ii)  Waves  with  large  enough  amplitudes  where  non-linear  effects. 


coupling  to  background  Brunt-Vaisala  Oscillation,  frequency  spreading,  induced 
instability  and  turbulence  become  important.  For  the  lower  altitude  waves  (m**dium- 
scale  T.I.D.'s),  where  many  of  the  sources  occur  near  the  ground  level,  the  increase 
in  velocity  amplitude  with  height  can  often  cause  the  Richardson's  number  to 
fall  below  critical  values,  thus  producing  instability  and  turbulence.  The  growth 
of  wave  amplitudes  will  eventually  be  damped  out  by  viscosity,  heat  conduction 
and  ion  drag  at  higher  altitudes.  The  wave  saturation  usually  occurs  above 
mesospheric  altitudes  which  leave  the  mesospheric  regions  with  the  largest  wave 
amplitude  and  hence  the  maximim  degree  of  turbulence. 

With  this  summary  of  gravity-wave  characteristics  we  see  immediately  that 
since  most  of  the  more  intense  infrared  emissions  come  from  lower  altitudes  (below 
100  km),  only  the  type  of  gravity  waves  which  produce  the  medium-scale  T.I.D.'s 
can  have  an  important  effect.  These  waves  have  been  investigated  by  Battaner 
&  Molina  (1980)  and  Hines  (1960).  The  results  may  be  summarized  as  follows: 

(i)  Gravity  Waves  vary  in  frequency  and  periods  in  the  range  between 
1  and  7  hours  with  the  most  often  occurring  periods  between  2  and  3.5  hours. 

(ii)  Seasonal  variations  show  a  Winter  and  Summer  maxima. 

(iii)  Horizontal  phase  velocities  vary  between  30  and  140  m  sec  * 
with  a  median  of  between  80-100  m  sec 

Based  on  the  above  summary  we  have  chosen  for  this  report  gravity  waves 
with  periods  of  two  hours  and  horizontal  phase  velocities  ranging  from  60  m  sec  ^ 
to  129  m  sec  Vail  of  which  correspond  to  the  most  frequently  occurring  low- 
altitude  gravity  waves. 

The  major  portion  of  this  report  deals  with  the  effect  of  linear  gravity 
waves  on  the  spatial  and  temporal  structure  of  minor  atmospheric  constituents 
which  produce  infrared  emissions  as  well  as  the  infrared  emission  profiles.  Since 
chemical  equilibrium  has  been  universally  assumed,  it  is  the  effect  of  gravity 


waves  on  the  initial  atmospheric  components  which  produce  the  excited  vibrational 
states  that  primarily  determine  the  spatial  and  temporal  behavior  of  the  infrared 
emissions  which  occur  as  a  result  of  these  vibrational  states  returning  to  the 
ground  states. 

We  have  found  two  possibilities  which  determine  the  basic  features  of  the 
spatial  structures  of  the  infrared  emission  profile: 

(1)  If  the  excited  vibrational  states  of  the  infrared  emitter  occur  locally 
as  a  direct  result  of  chemical  reaction  such  as  OH  emission  which  occur  mainly 
due  to  the  exchange  reaction  of  hydrogen  with  ozone  leaving  the  hydroxyl  molecule 

in  an  excited  vibrational  state  (Myrab<{>  et  al  (1983)),  the  spatial  infrared  profile  would 
show  wave-like  structures  just  as  strongly  as  the  number  density  profiles  of 
hydrogen  and  ozone. 

(2)  If  the  excited  vibrational  states  are  not  produced  locally  such  as 
the  case  for  chemical  reactions  but  are  mainly  produced  through  absorptions  of 
infrared  radiation  from  extended  sources  such  as  either  the  earth  considered 

as  a  blackbody  radiation  source  or  other  atmospheric  molecules  which  extend  over 
the  entire  atmosphere,  the  structures  would  be  modified  through  integrations 
over  the  extended  sources.  The  infrared  emitters  such  as  CO^,  NO,  0  and  H^O 
considered  in  this  report  largely  belong  to  this  class.  For  such  emitters  the 
infared  radiation  profile  would  be  less  sharply  structured  than  the  individual 
number  density  profiles. 

The  temporal  behaviour  would  in  general  be  far  less  strongly  influenced 
by  the  nature  of  the  source  and  the  excitation  mechanism.  Since  part  of  the 
present  report  is  in  collaboration  with  Visidyne  which  is  at  present  working 
on  the  temporal  radiation  profiles,  we  will  postpone  discussion  of  this  to  a 
later  report . 

The  general  characteristics  of  wave-induced  infrared  spatial  structures 


must  show  similar  characteristics  as  wave  structures.  Thus,  we  would  expect: 

(1)  Spatial  structures  relative  to  background  should  increase  with  altitude, 
since  such  structures  depend  generally  on  gravity-wave  velocity  amplitudes  which 
also  increase  with  altitude.  Irregular  structures  as  seen  in  the  CO^  radiance 
profile  by  Sharma  and  Nadile  (1980)  are  unlikely  to  be  due  to  a  gravity  wave 
because  the  structure  occurring  between  95  km  and  110  km  does  not  repeat  itself 
with  greater  severity  at  higher  altitudes.  However,  if  a  gravity  wave  has  high 
horizontal  phase  velocity  (and  hence  long  vertical  wave  length)  it  is  possble 

that  a  single  structure  will  appear  at  say  110  km  but  will  only  reappear  at  sufficiently 
high  altitudes  that  damping  would  serve  to  erase  it.  In  fact,  precisely  such 
a  structure  appears  near  110  km  in  the  CO^  profile  for  a  horizontal  phase  velocity 
of  129  m  sec  *  in  Fig.  8.  For  gravity  waves  with  shorter  vertical  wave  length 
(Fig.  13)  and  lower  horizontal  phase  velocity  (60  m  sec  ^)  we  see  that  the  structure 
increases  in  amplitude  further  up.  However,  the  damping  does  not  really  become 
excessive  until  several  hundred  kilometers  and  we  hence  regard  the  possibility 
that  the  structure  seen  by  Sharma  and  Nadile  as  being  caused  by  a  gravity  wave 
as  small. 

(2)  If  the  minor  atmospheric  constituent  has  a  layered  structure,  in  other 
words  it  is  not  in  hydrostatic  equilibrium,  the  oscillation  of  the  constitutent 
should  be  180°  out  of  phase  with  the  gravity  wave  where  its  ambient  concentration 
gradient  is  positive,  opposite  to  the  major  atmospheric  constituents  (bottom 
side  of  layer)  and  in  phase  with  the  gravity  wave  where  its  ambient  concentration 
gradient  is  negative,  the  same  as  the  major  constituents  (top-side  layer).  This 
result  has  been  shown  by  Hatfield  et  al  (1981)  and  has  a  long  history.  The  corresponding 
infrared  emission  should  show  a  similar  phase  change,  especially  if  the  excited 
vibrational  states  are  produced  locally  from  a  chemical  reaction. 

(3)  We  have  investigated  the  infrared  variance  with: 


(a)  Free  gravity  waves  with  vertical  wave  length  short  compared  with 
background  atmospheric  structure. 

(b)  Guided  gravity  waves  (Lamb  Mode)  with  vertical  wave  length  long 
compared  with  background  atmospheric  structure  e.g.  Fig.  2)  We  have 
found  that : 

(i)  Gravity  waves  with  short  vertical  wave  length  can  produce 
both  spatial  and  temporal  variations  in  the  number  density 
profile  and  hence  in  the  radiance  profile. 

(ii)  Gravity  wave  with  long  vertical  wave  length  such  as  the 
guided  Lamb  mode  can  only  produce  temporal  variations  in  the 
number  density  and  also  the  infrared  emission  profiles.  For  this 
case  the  radiance  profile  should  show  no  extraneous  spatial 
structures  (Fig.  2  and  3)  but  only  periodic  variations  as  a  function 
of  time  in  the  zenith  columnar  intensity  in  infrared  emissions. 

The  remainder  of  the  report  deals  with  gravity  waves  with  large  enough  amplitudes 
to  produce  instability  in  the  mesospheric  regions.  We  have  found  that  for  kz  X  >  1 
(where  kz  is  vertical  wave  vector  and  X  is  the  vertical  displacement  amplitude) 
which  typically  corresponds  to  X  ~  1.5  km  in  the  mesospheric  region  instability 
can  occur  and  give  rise  to  secondary  oscillations  with  periods  of  the  order 
of  Brunt  period.  In  very  rough  calculations  we  found  that  X  =  1.5  km  in  displacement 
would  require  a  density,  pressure  or  temperature  variation  of  20%.  All  the  experimental 

I 

evidence  in  airglow,  (Okuda  (1975)),  P.R.D.  (Manson  et  al  (1981))  and  temperature 
waves  (Myrabo  et  al  (1983))  seem  to  suggest  that  such  large  variations  will  produce 
secondary  oscillations. 

The  spatial  analyses  of  these  secondary  oscillations  at  an  altitude  close 
to  the  mesopause  and  with  periods  close  to  the  Brunt  period  cannot  be  done  with 


the  Hines'  dispersion  relation.  A  new  dispersion  relation  which  remains  valid 
in  an  inhomogeneous  medium  has  been  developed  for  this  report.  The  results  show 
that  the  Hines'  dispersion  relation  holds  well  for  all  but  the  frequencies  close 
to  either  the  Brunt  or  the  acoustic  cut-off  frequencies.  Close  to  such  frequencies 
the  vertical  wave  vector,  kz,  tends  to  zero  near  the  turning  point  where  heavy 
wave  reflection  takes  place  and  the  Hines'  relation  derived  for  a  homogeneous 
isothermal  medium  can  no  longer  be  valid.  It  is  also  precisely  within  the  range 
of  these  frequencies  that  all  the  observed  secondary  oscillations  appear  to  lie. 

Hence,  the  generalized  dispersion  relation  has  been  developed  to  analyze  secondary 
waves  arising  from  these  secondary  oscillations. 

The  detailed  technique  required  to  analyze  the  optical  emissions  data,  the 
temperature  wave  data  as  well  as  to  extract  from  these  data  the  means  to  determine 
atmospheric  parameters  such  as  turbulence  viscosity  and  energy  dissipation  rate 
will  be  relegated  to  a  future  report.  Experimental  procedures  for  measuring 

i 

long  term  infrared  optical  emissions  in  order  to  obtain  morphology  on  long-term 

I 

mesospheric  disturbance  will  also  be  part  of  this  future  program. 
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